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Abstract 

We exhibit explicit expressions, in terms of components, of discrimi- 
nants, determinants, characteristic polynomials and polynomial identities 
for matrices of higher rank. We define permutation tensors and in term of 
them we construct discriminants and the determinant as the discriminant 
of order d, where d is the dimension of the matrix. The characteristic 
polynomials and the Cayley-Hamilton theorem for higher rank matrices 
are obtained there from. 

1 Introduction 

A matrix A of rank r is an array of numbers Ai-^...i^ where the indices i run 
from 1 to d, the dimension of the matrix. 

For r = 2 we have the ordinary matrices a, that is, square arrays of numbers 
ttij. In this case it is possible to define a matrix addition and a matrix mul- 
tiphcation. Then, according to the Cayley-Hamilton theorem, only the first d 
powers of a are linearly independent. The traces of these d powers are a set of 
invariants which help to characterise the matrix a and several of its properties. 

For r > 2 the situation is different since, even when there is a matrix addi- 
tion, a natural matrix multiplication does not exist. The absence of a natural 
multiplication operation makes difficult to introduce concepts analogous to those 
which are standard in matrix calculus, namely, traces, invariants, polynomials, 
etc. Therefore, for matrices of higher rank there is not a result similar to the 
Cayley-Hamilton theorem. Furthermore, in principle, we have neither a way of 
determining the number of algebraically independent invariants (if any), nor a 
way to construct polynomial identities. 



Several results related to invariants and polynomial identities for matrices of 
higher rank are known [3 [9]. However, in general, explicit expressions in terms 
of the components Ai-^...i^ are lacking. 

The purpose of this article is to show that, in spite of the absence of a matrix 
multiplication, it is possible to define invariants, in fact a finite number d, that 
it is possible to define a determinant, and that there exists a polynomial identity 
similar to the statement of the Cayley-Hamilton theorem. 

In section 2 we start with a reminder of some standard results for ordinary 
matrices. In this case, invariants can be obtained, for example, as traces of 
the powers of the matrix a. Other family of invariants are the discriminants. 
Discriminants are a more convenient set of invariants than traces since only the 
first d of them are non trivial, while the rest are identically zero. Particularly 
interesting is the discriminant of order d which corresponds to the determinant 
of a. The Cayley-Hamilton theorem is formulated as the condition for the 
vanishing of a certain polynomial relation among discriminants and powers of 
the matrix a. Next, we introduce an index notation similar to that of tensor 
calculus and all matrices are then treated as tensors. Discriminants can be 
defined in terms of permutation products with a metric tensor g. Since, this 
second definition makes no reference to traces, this definition of discriminants 
is best suited than traces for the generalization to higher rank matrices. 

In section 3 we consider higher rank matrices. Since there is not a concept 
similar to that of a multiplication operation, the construction of invariants must 
be done by generalizing definitions which are independent of the concept of 
multiplication. For instance, discriminants can be generalized to higher rank 
matrices if we adopt the definition based on permutation products. The even 
rank and the odd rank cases must be studied separately. For even rank matrices 
the fourth-rank matrices are the prototype. By using permutation products 
we construct the corresponding discriminants and define the determinant as 
the discriminant of order d. Then we show that a certain polynomial relation 
among discriminants and products of the matrix A vanishes identically, which 
is a statement similar to the Cayley-Hamilton theorem. For odd rank matrices, 
a naive generalization of the results above leads to useless relations. Instead, 
we show that in this case it is necessary to introduce an even rank matrix as 
the direct product of the original odd rank matrix and construct discriminants 
and polynomial identities as for the even rank case. 

Section 4 is dedicated to the conclusions. For the sake of simplicity in the 
present work we restrict our considerations to completely symmetric matrices 
and tensors. Generalizations to other situations (other ranks and more general 
symmetries) are easy to implement, but we refrain to exhibit them for they do 
not add any new understanding to the problem. 

Some preliminary and imcomplete results, similar to the ones reported here, 
can be found in 20J. In that work our emphasis was on the construction of 
invariants using a graphical algorithm based on semi-magic squares. However, 
the rigorous mathematical justification for that algorithm are the invariants 
constructed with matricial properties alone as exhibited here. 
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2 Preliminaries 



Due to the nature of our approach, we must regretfully bore the reader by 
exhibiting some standard and well known results that are needed to show how 
the concepts of invariant, discriminant, determinant, characteristic polynomials, 
polynomial identities, and other concepts of matrix calculus, correctly generalize 
to matrices of a higher rank. Even when these results can be found in standard 
references, see for example [5], we include them here because we need them 
written in very particular forms (which are not to be found in a unique reference) 
which make easier, even conceptually, the passage to matrices of higher rank. 



2.1 Basic definitions and results 

A matrix a is a square array of numbers a^ , with i.j — 1, • ■ • , d, where d is 
the dimension of the matrix. Let b be a second matrix with components bij. 
Then, the matrix multiplication is defined by the product matrix c = a ■ b with 
components 



= ark bkj ■ (1) 



fc=i 



This operation is the Cartesian product among rows and columns. The unit 
element I for the matrix multiplication, a • I = I • a = a, is the matrix with 
components hj given by 

J . ^fl, fori=j; 
[ , otherwise. 

The determinant of a, det(a), given by the Leibniz formula 

det(a) = sign(7r) Y[ a„(j) j . (3) 

If the determinant ([3]) is different from zero, then there exist an inverse matrix 
a~^ satisfying ■ a — a. ■ = I. In terms of components we have 



(a ^)ik akj = Y "^'t ^ ■ ('^) 

k=l k=l 

Then the inverse matrix is given by 

1 _ adj(a) 

where adj(a) is the adjoint matrix. 

The product of a matrix a with itself, that is a^, is the matrix with compo- 
nents 
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d 

(a^)y = ^ aik akj ■ (6) 

fe=i 

Powers of a of an order s, a*, are the matrices with coraponents 

d d 

= X! "' «»fci • • • "fc,_ij ■ (7) 

fcl = l fes-l = l ^ ' 



s — 1 times 

By definition a'' = I and a^ = a. The trace of a matrix a is given by 

d 

tr(a) = ^ a,, . (8) 

i=l 

The trace of a* is given by 

d 

tr(a^) = ^(a^k. (9) 

1=1 

Furthermore, tr(a") = d. 

Let us now consider the eigenvalue problem 

[a-AI]v = 0. (10) 
The condition for the existence of a solution is 

det[a-AI] = 0. (11) 
From here it follows that the eigenvalues are related to the traces by 

d 

tr(a«) = 5] A^ (12) 

k=l 

An invariant is a quantity which does not change under similarity transfor- 
mations. Invariants are important because they help to characterise a matrix a 
and several of its properties. Examples of invariants are the traces Q. 

Discriminants are another family of invariants which are constructed as fol- 
lows. Let us consider a vector A = Ai, A2, • • ■, with an infinite number of com- 
ponents. Then, the elementary symmetric products P are defined by 

Ps{\)= ^k,---Xk^- (13) 

The power sums Q are defined by 

00 

Q.(A) = 5] A^ (14) 

fe=i 
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Elementary symmetric products and power sums are related by 



Pq 


= 1, 


Pi 


= Qi, 


P2 




Pz 


- ^[^? 


Pa 




P5 





- [Qt ~ 6QIQ2 + SQ1Q3 + 3Ql - 6Q4 , 

- 10 Ql Q2 + 15 Qi Ql + 20 Ql Q3 
2OQ2Q3-3OQ1Q4 + 24Q5] , 



- [Ql - 15 Qt Q2 + 45 Ql Ql - 15 Q: 



+40 Ql Q3 - 120 Ql Q2 Q3 + 40 Ql 

-90 Ql Qi + 90 Q2 Q4 + 144 Qi Q5 - 120 Qg] 



(15) 



These are the Newton relations. 

If only the first d components of A are different from zero, then Ps = for 
s > d. 

Therefore, traces are in correspondence with elementary symmetric products. 
In order to complete the correspondence, and following relations p^ . the dis- 
criminants of a matrix a are defined by 



C6(a) 



Co (a) 


= 1, 










ci(a) 


= (a), 








C2(a) 


1 - 

2 . 






7 




C3(a) 


1 
3! 




-3(a)(a^) + 2(a3); 




C4(a) 


1 
4! 




- 6 (a)' (a^) + 8 (a) (a^) + 3 (a^)' - 6 (a^)" 


C5(a) 


1 
5! 


;(a)' 


- 10 (a)' (a^) + 15 (a) (a^)' + 20 (a)' (a^) 



-20(a2)(a3)-30(a)(a'*)+24(a^)] , 



1 
6! 



(a)° - 15 (a)" (a^) + 45 (a)' (a^)' - 15 (a^) 

+40 (a)' (a3) - 120 (a) (a^) (a^) + 40 (a^)' 
-90 (a)' (a^) + 90 (a^) (a^) + 144 (a) (a^) 
-120 (a^)] , 



(16) 
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etc., where (•) = tr(-). Discriminants of a square matrix are not so widely 
known objects. They appear in the htcrature under different names or indirectly 
in several contexts. 

The discriminants (|16p satisfy the remarkable recurrence relation 



^^^^"^ a-c.(a)I.-^^. (17) 



da. dsL 
where {dcs (a) / da) is the matrix with components 

dcs{a)\ dcs{a) 



da J duij 
We can state the Cayley-Hamilton theorem as: 

Theorem. (Cayley-Hamilton) A d-dimensional matrix a satisfies 

dcd{a) 



da 



(18) 



a-Cd(a)I = 0. (19) 



This result follows from (fT7)l . for s = d, reminding that Cd+i(a) = 0. For the 
first values of d the explicit statement of the Cayley-Hamilton theorem is 



a — ci (a) I = , 
a^ -ci(a)a + C2(a)I = 0, 
a^ --ci(a)a2 +C2(a)a-C3(a)I = 0, 
a''-ci(a)a3+c2(a)a2-C3(a)a + C4(a)I ee 0, (20) 

etc. 

According to the Cayley-Hamilton theorem, only the first d powers of a 
are linearly independent. Therefore, only d of the traces ([9]) are algebraically 
independent. The ideal situation would be therefore to have a family of invari- 
ants such that only d of them are non-trivial and algebraically independent. 
Discriminants have this property. For a d dimensional matrix, only the first d 
discriminants are non-trivial, while the discriminants of an order higher than 
d are identically zero, Cs{a) = 0, for s > d. This result is equivalent to the 
Cayley-Hamilton theorem. This reformulation of the Cayley-Hamilton theo- 
rem is possible because we can establish an equivalence between traces, power 
sums, elementary symmetric products and then Cs = 0. Therefore, we consider 
discriminants as our fundamental family of invariants. 

If Crf(a) 7^ 0, then from ([19]) follows that there exist an inverse matrix a^^ 
which is given by 

a- -4^^. (21) 
Cd{a) da 

In terms of components this inverse matrix is given by 
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Cd{a) da,j 

We have two algorithms to compute the inverse of a matrix a: 

Al The first one is based on ([4]). In that case we have n? unknowns (a~^)ij 
and n? equations. The condition to have a solution is det(a) ^ 0. 

A2 The second algorithm is based on the discriminant Cd(a). If Cc((a) ^ 0, 
then the inverse matrix is given by (1211) . 

If Al, • • ■ , Ad are the eigenvalues of the matrix a, then Cd{a) Ai ■ • • A^, which 
is the definition of the determinant of a. Therefore, Cd(a) = det(a). This result 
explains why Cc;(a) and det(a), even when constructed with different algorithms 
lead to the same inverse matrix, namely ^ and ([2T|) . For simplicity we denote 
a = det(a) = Cd(a). Then, the equation (PT|) is rewritten as 

a- = -|^- (23) 
a oa 

For ordinary matrices, a, the two previous algorithms are equivalent and yield 
the same inverse matrix a~^. However, for higher rank matrices only the second 
algorithm, based on (pij) . admits a generalization. 



2.2 Matrices and tensors 

For the purpose of generalization to higher rank matrices it is convenient to 
represent matrices by means of tensors. We adopt the index notation of tensor 
calculus. In this case it is necessary to distinguish between covariant and con- 
travariant indices. According to this scheme, a matrix a becomes represented 
by a second-rank covariant tensor a with components . Also we will need to 
consider contravariant tensors g with components g^^ . Furthermore, we adopt 
the summation convention according to which an index is summed up over its 
range if it appears once as a covariant index and once as a contravariant index; 
for example, aikg^^ means X]fc=i o.ikg^'^ ■ 

Let a be a covariant tensor with components . The unit element e for the 
tensor multiplication is a tensor with components e'j given by 

^ *[ , otherwise. ^2'^) 

The inverse a~^ is a tensor with components (a~^)*-' such that 

{a-^r'' a,k ^ 5] ■ (25) 

Since the inverse tensor is a contravariant tensor the notation [a^^Y^ becomes 
redundant; therefore, we simply write a'-' for the components of the inverse 
tensor. For example, in term of components, equation (|23|) becomes 
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1 da 

- 7^ , 26 

a oa.ij 



while equation (|25|) becomes 



a^^'ajk^S). (27) 

In order to define a multiplication operation in a tensor language it is neces- 
sary to consider a metric tensor g with components gij. If g ~ det(g) 7^ 0, then 
there exist an inverse with components g'-' . Then, if a and b are covariant 
tensors with components Uij and 5^ , the multiplication operation is defined by 
the product tensor c with components 



The trace is defined as 



a,k g"" bij . (28) 



trg(a) = .9*^ fly . (29) 



The definition of the trace in terms of the same metric tensor g used for the 
multiplication is necessary in order to preserve the property of ciclicity of the 
trace. 

In order to agree with ([T]), it would be necessary to choose 

^ "1 , otherwise. ^^^^ 

However, this definition has not an invariant tensorial meaning and therefore it 
is better to keep working with a generic metric tensor g^^ with components (7'^ 
unrelated to (l30l). 



2.3 Permutation products and discriminants 

For higher rank matrices neither a natural multiplication operation nor a natural 
extension of the concept of a trace do exist. Therefore, we need to elaborate a 
definition of discriminants without making reference to multiplications or traces, 
and, in such a way, that the resulting definition reduces to the usual definition 
for matrices of second rank and is easily generalizable to matrices of higher 
rank. To this purpose, we define permutation products in terms of g. As a 
motivation and justification to do so we start by considering the definition in 
tensor calculus of the determinant of a tensor a. 
The Levi-Civita symbol is defined by 

{1 if ii • ■ • id is an even permutation of 1 • ■ • d; 
— 1 if ii • ■ • id is an odd permutation of 1 • • • d; (31) 
if ii ••• id is not a permutation of 1 ■ • • d. 

Then, the determinant of a tensor a with components Uij is given by 
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det(a) = le'--'''e^--^"-a,,,, •••a,,,,. (32) 
Substituting ([5^ in we obtain the explicit expression 

We can verify that a'-' so defined satisfies (P7|. It can also be verified that 
equation p3p is equation ([5]) written in terms of components. 

Therefore, for second rank covariant tensors the discriminants can be con- 
structed with the use of g for products and traces, as in and (P^ . and 
the definitions for discriminants in (jl6p . However, for the determinant, which 
corresponds to the discriminant Cd(a), we have (|32p which involves a alone and 
that makes no reference to g. In order to avoid, and explain, this duplicity we 
develope a unified scheme to construct discriminants and which is well suited 
for the generalization to higher rank matrices. 

Let us consider a metric tensor g with componentes gij and define the per- 
mutation tensors q by 



9s ^^'^ ~ s!(d - s)! 9is+ijs+i yidjd- 

The tensors q are non trivial only for s < d, and q^ = for s > d. We define 
the discriminants cf (a) for a tensor a by 

cf(a) = g:^^--»^»(g)a,,,, •••a.,,,. (35) 
Then, the discriminants cf (a) are given by 



^(^\ — ^ ,»i---'3»3 + i---'d ,.7i---.?3.?3 + i---.7d „ n n n 

(36) 

For second-rank tensors, the case under consideration, the permutation ten- 
sors (p4| can be rewritten as 

ql^n■■■^.Js^^)^]_g\l^^n (37) 
31 

where |[- • •]! denotes complete antisymmetry with respect to the indices j's or, 
equivalently, with respect to the indices i's. For the first values of s the tensors 
q are given by 
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^m23...0. (g) ^ 1 [^.Ul g^.3. gr.J. _ ^gr.n g^.n g^.3. 



+ g..j3 g^sn + ^nj3 g.3j.)] ^ (38) 

etc. When we restrict g to ([30|) . (|36| reduces to the relations (fT6|l with (•) = 
trg(-); therefore this is a valid generalization of the discriminants. 
We can verify that 



9(gcf(a)) ^ g(gcf^i(a)) 
and from here it follows that 



(39) 



^+<f(a)g-..%<5>, (40) 

This is a recurrence relation analogous to ^T7\ . 

We can now reformulate the Cayley-Hamilton theorem as: 

Theorem. (Cayley-Hamilton) A d-dimensional tensor a satisfies 



MM 

9g 



c^(a)g-i^O. (41) 



This results follows from (I40p . for s — d, reminding that c^_^^{a) = 0. 
There are two particularly interesting instances of relation (|4ip : 

Bl The first case appears if we insist on g~^ as given in ([501) attempt 
to reproduce the standard results. Then, the relations in (|1T|) reduce to 
the relations of matrix calculus given by (|20p . 

B2 The second case appears if we want to have expressions concomitant of a 
alone. In this case the only possible choice is g = a. However, all relations 
in (|4T|) collapse to useless identities. 



The generalization of the first case to higher rank matrices is excluded since it 
is not possible to define a metric tensor g similar to an identity or unit matrix 
in an invariant way. On the other hand, the generalization of the second case 
is possible and is the one which gives rise to discriminants and the Cayley- 
Hamilton theorem for higher-rank matrices. 

The componentes of are completely antisymmetric in the d indices i and 
the d indices j. Therefore, they must be the product of Levi-Civita symbols in 
those indices 

^n---idji---jc!(g-j _ ]_ g\[iiji . . . gidjd]\ _ 1 i_ gil gii-'-jd _ (42) 
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Then, we obtain 




a 



9 



(43) 



Let us denote ag = c^(a) 



a/9 



If a 7^ and g ^ 0, we obtain 



a 



1 



1 da 
a da. 



1 dag 
ag 9a 



(44) 



This relation explains why the determinant, as defined in p2p (which does not 
involves g), gives the same inverse tensor as c^(a). 

3 Higher rank matrices 

The concept of a higher rank matrix, and the corresponding determinant, was 
introduced by Cayley in 1845 [5] and it was later developed by Schlafli [T5] 
and by Pascal [TO]- More recently, matrices of higher rank have been studied 
in [71 [SI [HI [131 [m [2S] . Particularly interesting are |S] and [S] where a general 
account on the subject, with many generalizations and applications, can be 
found. The interested reader may consult these references for further detail. 

For higher-rank matrices there is not a natural multiplication operation in 
the sense that the product of two higher rank matrices be a matrix of the same 
rank and covariance. Therefore, the construction of discriminants must be done 
by generalizing a definition of discriminants which do not involve any multi- 
plication operation. The definition of discriminants in terms of permutation 
products satisfies this requirement, but even rank and odd rank cases must be 
considered separately. For the even rank case we take the fourth-rank case as a 
prototype. Then, we construct discriminants, the determinant, and show that a 
polynomial relation among discriminants and products of the higher rank ma- 
trix vanishes identically, which is a statement similar to the Cayley-Hamilton 
theorem. This algorithm can be easily extended to matrices of any arbitrary 
even rank. For the odd rank case we consider third-rank matrices as a proto- 
type and show that in this case it is necessary to introduce an even rank matrix 
constructed as the direct product of the odd rank matrix. Then, discriminants, 
the determinant and the Cayley-Hamilton theorem follow as for the even rank 
case. 

Let us remind once again that higher-rank matrices will be represented by 
higher-rank tensors. Therefore, we will talk, mostly, of tensors rather than 
matrices. 

3.1 The even rank case 

We consider fourth-rank matrices as a representative of even-rank matrices; all 
other even-rank cases can be dealt with in a similar way. 
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The simplest extensions of the concept of determinant to matrices of higher 
rank is the one due to Cayley, which is a direct gencrahsation of the Leibniz 
formula; namely 

detc(A) = sign(7r2 • • • tt^) \Yl A (45) 

7r2---7r,.e5'd \i=l / 

Then, the determinant of a fourth-rank tensor A with components Aij^i is 
defined as a direct generahzation of ([5^ . namely, 

det(A) = 1 e^--^- ■ ■ ■ e'--'^ ■ ■ ■ A,,,,^,,, . (46) 

As before we simplify the notation by writing A ~ det(A). In analogy with (|23p 
we define 



In term of components 



Then 



(48) 



ijkl 



Aijkl _ ^ _J_ Jii---i(d-l) . . . Jll---l{d-l) 4 4 . 

^ ~ (d— 1)! A ^njifeih ^i(d-i}J{d-i)k(d-i)l(d-i} ■ 

(49) 

This higher rank tensor satisfies, 

A'"'"'''^ A,kM = 5] , (50) 

which is a relation similar to (|27p . The definitions in and were used in 
previous works [121 [U [H] concerning the application of fourth-rank geometry 
in the formulation of an alternative theory for the gravitational field. 



As an example let us consider the simple case d = 2. The determinant (|46p 
is given by 

A — Aoooo Ann - 4 Aqooi Aqui + 3 Aqqh . (51) 
The components of the inverse matrix A"-'^ are given by 

40000 _ 1 4 
A°°"i = -^Aoni, 

A™" = ^Aoon. (52) 



12 



and similar expressions for the other components. In order to verify the vahdity 
of equation (jSO]) let us consider the cases (00) and (01). We obtain 



A°'^'Ao„fc = 1, 

A°'^'AH,fc = 0, (53) 
and similar relations for the other indices. 

In a way similar to (j34p we define the permutation tensors Q by 

^ ' ' s\{d-s)lG 

Gis+ijs + iks + ih + i ■ ' ' ^idjdkdld 1 (54) 

where G = det(G). The tensors Q are non trivial only for s < d, and Qs = 
for s > d. Then, we define the discriminants C^(A) by 

(A) = Ql^^^'^'--^^^'^'^{G)A,,,kj, ■ ■ ■ A^j^kj^ . (55) 

Explicitly 



r<'^(A \ — ^ ^ Jl-'-isiB + f-id Jl---jsjs+l---jd c^l---ksks + i---kd J-i-'-lsls+l-'-ld 

^^^^>- s\{d-s)\G ' ' 

X ^iijifci/i • • • Ai^j^k^i^ Gi,^_ij3_|_ife^+i/^+i • • • Gijj^kdid ■ (56) 

Then, in a way similar to (|40p . it is possible to verify that the discriminants 
Cf (A) satisfy the recurrence relation 

dCfjA) ggf+i(A) 

+ (A)G = . (57) 

Then, we have 

Theorem. (Cayley-Hamilton) A d- dimensional tensor A of fourth-rank sat- 
isfies 

^^ + Cf{A)G-'^Q. (58) 

The components of Qd are completely symmetric in indices i, j, fc, and I; 
therefore, they must be the product of Levi-Civita symbols. We then have 

Qiljlkilf-idjdkdld I^Q^-^ _ }_ }_ ^il---id ...^If-ld ^gg^ 

G d\ 
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which is the result similar to (|^^ . Therefore 



Cf{A) = ^. (60) 
Let us denote Aq = Cf{A) = A/G. II A^Q and G 7^ 0, we obtain 

A-i-iM-^M£ (61) 



which is the relation similar to (j44| . 

There is one particularly interesting instance of the relations ([58|) . If we 
want to have expressions concomitant of A alone, then the only possible choice 
is G = A. To this purpose it is useful to observe that if s = d, and only in 
this case, the permutation tensors (j54p can be written in a way similar to (|37p . 
namely 

Qiljlkilf-idjdkdld _ }_Q\[iljlkili _ ^ Qidjdkdld]\ (g2^ 

As an example of the results above let us consider the case d = 2. For the 
tensor Q2, we obtain 

Qhjikilii2j2k2l2 ^Q,-^ _ }_ ^Qiijikih Qi2j2k2l2 ^Q'hjikih Qi232k2li 

_^Qiljlk2ll Qi2j2kll2 _|_ Qilj2klll Qi2jlk2l2 
_^Qi2jlklll Qilj2k2l2^ _^ ^Qiljlk2l2 Qi2j2kili 
_^Qilj2kll2 Qi2jlk2ll _j_ Qi2jlkll2 Qilj2k2ll^'j (63) 

The determinant is given by 

(A) = i [(G'^"'=' A,,ki)^ ~ 4 G*^''^' A.kim G'^^P" A„p,, 

+3G'^'='A«„„G™Ap,„], (64) 
Finally, the corresponding polynomial identity is given by 

f^Qirmpq Amnpq) Aijkl — 4 A(j|„j„p Qmnpq j^^^^.^^^ 

+3A(ij|„„G™"'"' Apqifej) - G^(A)Gij7ci = 0. (65) 
If we choose G = A the relation (ISS)) reduces to 

■^{ij\mn A"^^P'^ ^jjqlfci) ~ ^ {A"^^^'^ Apgrs A""**" Atumn) Aijkl = , (66) 

where (•) means that the enclosed indices are symmetrised and | • | means that 
the enclosed indices are excluded from the symmetrisation. 
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3.2 The odd rank case 



We consider third-rank matrices as a representative of odd-rank matrices; all 
other odd-rank cases can be dealt with in a similar way. 

Let s be a third-rank tensor with components Sy j,. A naive definition of the 
determinant would be the natural generalization of (|32p and (|46p . namely 

det(s) = 1 6^--^^ . . . e'^^-'^- s.,,,,, ■ • ■ s,,,,,, . (67) 
However, it is possible to verify that 

e''-'' ■•■^'=-■■'=^5,,,,,, ■■■ s,,,,k,^0. (68) 

This result is due to the odd number of e's in (|S7|) which add all contributions 
to zero. 

In order to obtain some indication as to the correct way to define a deter- 
minant for odd-rank tensors, let us consider the simple case of a completely 
symmetric third-rank tensor s with components Sijk in dimension d. Then, 
let us look for an inverse tensor s^^ with components s'-'*^ such that a relation 
similar to ([TT]) and ([50)1 holds, namely 

s'"' sjki = 5] ■ (69) 

The number of unknowns in (|69p is d{d + l){d + 2)/6 while the number of 
equations is c?^. This algebraic system of equations is underdetermined, except 
for d = 2. In this last case the solution is given by 



s°°° = ^ (soooSiii + 2sqii - 3sooi Son Sill) , 
s°°^ = ^ (2 Sill Sqoi - sooo Son Sill - sooi Soil) , (70) 
and similar relations for the other components, where 



Q 322/\ 

s = Sqoo Sill - 6 sooo sooi son sm + 4 sooo Son + 4 sm Sgoi - 3 Sqoi Squ . (71) 
Let us observe that 



s°"° 



1 (9s2 Id 



s 



2s2 dsQoa s dsooo ' 
3 2s^ 9sooi 3 s 9sooi ' 



leading to 



s 



ijk 



1 ds^ 1 ds 

2s2 dsi-jk s dsijk 



(73) 
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Therefore, the role of the determinant, in a way similar as appears in (j26p and 
(|48l), is played by s. 

For any matrix S the product of e's and S's is a meaningful quantity only 
for an expression in which S is an even rank tensor. Then must be related to 
the determinant of some higher even rank tensor S, and this determinant must 
be a quadratic expression (since d = 2) in this higher rank tensor S. Since 
is a quartic expression in s we have that S must be quadratic in s. A solution 
satisfying the requirements above is given by 

For d = 2, we obtain 



'S'oOOOOO — Sgoo I 

5*000001 = sooo sooi , 

5ooooii = - (2 sooo Son + 3 Sqoi) , 

5oooiii = ^ (sooo Sill + 9sooi Son) , (75) 

and similar relations for the other components. 

The determinant of a sixth-rank tensor S with components Si-^i^i^i^i^i^ is 
defined through a direct extension of the definition in P5|) . We obtain 



5 — — e*^ • • ■ e"^ "'^ 5ij...„i ■ • ■ 5i^...„^ . (76) 



1 

In the the two-dimensional case this determinant is given by 



5* — 5oooooo 5*111111 — 6 5oooooi 5*111110 + 15 5ooooii 5iiiioo — 10 Sqqq^h . (77) 

If we replace (|75p in (j77[) we obtain, up to an irrelevant multiplicative constant, 
the expression ((7T|) . 

Therefore, for odd rank matrices the recipe is to construct an even rank 
matrix as the direct product of the original matrix. Then, the construction of 
invariants proceeds as for the even rank case. 



4 Concluding remarks 

We have developed an algorithm to construct invariants for higher rank matrices. 
We constructed determinants and exhibit an extension of the Cayley-Hamilton 
theorem to higher rank matrices. 

Among the possible applications of the present work we have the quantum 
mechanics of entangled states. Indeed, in order to obtain a measurement of en- 
tanglement one needs to construct invariants associated to higher rank matrices; 
see [2 H H El [H m m m m HI] for several interesting results. 
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Higher rank tensors, which look similar to higher rank matrices, appear in 
several contexts such as in Finsler geometry [2j [17] and in fourth-rank gravity 
[ini mi [21] • The results presented here are a first step for the construction of 
differential invariants for higher-rank tensors. 
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